# -*- coding:utf-8 -*-
"""
Created on 2009-10-28

@author: Qiu wenfeng
"""
import scipy

def curvature_central(nx, ny):
    junk, nxx = scipy.gradient(nx)
    nyy, junk = scipy.gradient(ny)
    K = nxx+nyy
    return K


